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ABSTRACT 

The main experimental results on jet production at HERA are reviewed. A 
study of jet shapes shows that the internal structure of jets is well understood. 
The potential to accurately determine the strong coupling constant using jet 
rates is discussed. First attempts to extract the gluon density in next-to-leading 
order are presented. Although fixed order perturbative QCD is able to describe 
the shape of most variables associated to the hard subprocess leading to jets, it 
does not reproduce the absolute dijet rate over a large phase space region. 

1. Introduction 

HERA colliding 27.5 GeV positrons on 820 GeV proton offers an ideal testing 
ground for perturbative QCD in deep-inelastic scattering (DIS). 

Hard processes can be investigated as function of a variable scale, such that the 
transition from one phase space region to another or the interplay between two hard 
scales, like e.g. the photon virtuality or the transverse energy of a parton, can be 
studied. The center of mass energy of i/i ~ 300 GeV leads to a large phase space 
for hadron production and to the possibility to observe clean 'jets' in the hadronic 
final state. 

Jets are event properties relating the unobservable quarks and gluons to the mea- 
surable hadronic final state. Quantitative studies of jets require an exact definition of 
how to combine the spray of hadrons observed in the detector to jets. The definition 
should be suitable for the experimental analysis and for the theoretical calculation. 
Moreover, it should lead to infrared and collinear safe results, which do not change, 
if a low energetic particle is added or if a particle is split into two. Jets are usu- 
ally defined either by cone type algorithms maximizing the transverse energy flowing 
through a cone in the [t], 0) plane^or by algorithms successively recombining clos- 
est particles. Clustering algorithms iteratively merge pairs of particles until only a 
few well separated objects are left. The decision to assign a given object to a jet is 
based on a distance measure dij and a resolution parameter y^ut = dij /scale. Table |1] 

" The pseudo-rapidity is rj = — Intan (6'/2), where 9 is defined with respect to the incident proton. 
The proton defines the negative z axis, (j) is the azimuthal angle. 
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Fig. 1. Summary of the main characteristics of popular jet clustering algorithms. For the Kt 
algorithm the remnant is considered as a particle p with infinite momentum. In the JADE and 
W algorithm a pseudo-particle with longitudinal components calculated from the hadrons in the 
final state is introduced to account for the remnant. The combination procedure defines how the 
four-momentum pi — {Ei,pi) of two objects i and j are combined to an object k. 



summarizes choices for commonly used clustering algorithms. Also shown are various 
procedures for the addition of particle 4-momenta to obtain the jet four-momentum 
(recombination schemes). 

Traditionally cone algorithms have been applied in pp collisions where special care 
has to be taken to separate particles belonging to the beam remnant from particles 
associated with the hard subprocess. Cone algorithms are invariant under longitudinal 
boosts. In e~^e~ collisions, where the initial state is free of color charges and where 
every particle can be assigned to jets, clustering algorithms have been used. In DIS 
or pp collisions they can also be used if an additional particle to mimic the hadron 
remnant is put in by hand. 

At HERA, events with two jets in the central part of the detector plus the jet 
associated with the proton remnant (2-1-1 jet events) can be produced in a quark 
(g7 qg) or a gluon {gj qq) initiated hard subprocess (see Fig. 0). The 2 -|- 1 jet 
cross section is given byu: 

d'^a2+i , f dC 



dxdQ"^ "s(/^r) J —{Cqg{x/i,Zq,X,Q )-q{^,fIf)+Cqq{x/^,Zq,X,Q )-^(^,/i^)) (1) 

where as is the strong coupling constant, fir the renormalisation and /i/ the factori- 
sation scale. X and are the usual variables to inclusively describe deep-inelastic 
scattering (see Fig. ^ for definition). The coefficient functions Cqq and Cqg can be 
calculated in perturbative QCD. In addition to x and they depend^on the two 
variables ^ and Zq characterizing the short distance subprocess which is responsible 
for the dijet system. The fraction of the longitudinal proton momentum carried by 
the incoming parton ^, can be calculated in leading order as (LO) from the invariant 
mass of the jet system 

^ = x{l + s/Q') (2) 

''The dependence of Cqg and Cqq on the azimuthal angle between the lepton and the parton 
scattering plane is not given in the formula. 
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DIS variables: 
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x = Q^ /{2P ■ q) 
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Fig. 2. Feynman diagrams for the production of 2 + 1 jet events to first order of as in ep -collisions. 
q (P) denote the four- momentum of the photon (proton), ji (^2) is the four- momentum of one of 
the jets associated to the hard subprocess. 



The variable Zq is given by: 

P-Ji 



P-q 



cosri 



^j(l-cos^,) 



(3) 



where q (P) is the 4-momentum of the photon (proton) and 6* is the polar angle*^ 
of one of the jets with 4-momentum ji in the photon-parton center of mass frame. 
Ej and 6j are the energy and polar angle of the jets in the laboratory frame. The 
parton distribution functions in the proton g(^, Q^) and g{^, Q"^) absorb the collinear 
and infrared singularities which occur in the calculations of Cgq and Cgg, and they 
have to be extracted from experiments. Once measured at some scale Qq, they can 
be perturbatively evolved to any other scale, e.g. Q^, using the DGLAP equationJl. 
By convoluting them with the appropriate coefficient functions, they can be used to 
calculate the cross section in any other process. In this sense they are universal. 

In leading order of as, the coefficient functions diverge like Cgq ~ l/{zq{l — Zg)) 
and Cgg ~ — — Zg)) and have to be regulated by a cut-off. In Fig. | is 

shown how the cross section d?a / dxpdzg (where Xp = x / ^) increases towards low Zq 
and Xp and towards Zq = 1. It is clear that in the regions where the LO cross section 
diverges, large contributions from higher orders can be expected. 

The dijet cross section including higher order parton emissions are difficult to 
calculate using exact expressions for the coefficient functions. Since the phase space 
integrals cannot be solved analytically, numerical methods have to be applied. Re- 
cently several Monte Carlo integration programs computing the cross section to next- 
to-leading order (NLO) became available. Early attempts like projetI or DISJETQ 
used a semi-analytical approach which only allowed jets to be defined by the W al- 
gorithm with the JADE scheme^ and were restricted to a small set of calculable 
distributions. Moreover, in these programs approximations were made which turned 
out not to be valid over the full phase spaceS. The programs MEPJEtI, based on 
a phase space slicing methodi, and DISENli^, based on a method to subtract the 



■^The polar angle is defined with respect to the incoming parton direction. 
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Fig. 3. a) Differential cross section cPa/dxpdzq in first order of for x = 10~^ and ~ fO GeV^. 
The figure is taken from ref.H. b) Relative change in percept of the dijet rates when adding to the 
leading order cross section higher orders based on a leading logarithm parton shower model in 
the LEPTO Monte Carlo. The lines represent typical selection cuts on the scattered electron. In 
addition, a cut on Zq > 0.1 has been applied. The jet rate is based on the JADE algorithm for 
Vcut = 0.02. The figure is taken from rein. 



divergences from individual contributionJHI and on an analytical treatment of specific 
U divergent termJll, are flexible NLO Monte Carlo programs which allow any jet defi- 
nition scheme and arbitrary experimental cuts to be analyzed. For > 40 GeV^ 
and for jets defined by the kr algorithmlll the agreement of the 2 + 1 and 3 + 1 jet 
cross sections between the two programs is found to be on the 3% levelli. 

Higher order parton emissions and models for the transition of partons into 
hadrons are only included in Monte Carlo programs incorporating coefficient func- 
tions to leading order such as ARIADNEji, LEPTO0, HERWIG0, PYTHlAlli and 
RAPGApi. They model multi-parton emission by either parton showersil (PS) 
based on the DGLAP equations or by a chain of independently radiating colour 
dipolesi (CDM). HERWIG uses the cluster hadronisation modelH. In all other 
programs hadronisation is based on the LUND string mo deii as implemented in 

jetsetE. 

The relative change in percent of the dijet rates defined with the JADE algorithm 
for ycut = 0.02 when adding to the LO cross section higher orders based on a leading 
logarithm parton shower model (MEPS) is shown in Fig. ^ in the x -Q^ plane applying 
typical analysis cuts on the scattered electron. While for Q"^ > 100 GeV^, the 
relative change only exceeds 20% in the low x region, parton showers change the 
dijet rate by a large factor for 10 < < 100 GeV^ and 1 ■ 10"^ < x <2 - 10^^. 
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Fig. 4. a) Rate of 2 + 1 jet events as function of the invariant mass of the two jet system 
and the invariant mass of the hadronic final state W. Indicated are the phase space restrictions 
introduced when using the cone or the JADE algorithm, b) Fraction of gluon and quark induced 
processes leading to 2 + 1 jet events with Et > 5 GeV and — 1 < 77 < 2 defined by a Kt algorithm in 
the Breit frame as function of . The cross section has been calculated in NLO in the kinematic 
range a; > 1 • 10^^ and 0.15 < y < 0.65. 

The influence of parton emissions beyond leading order does also depend on the jet 
variables Zq and ^. Since different jet algorithms cover different phase space regions, it 
is always important to specify which jet definition scheme was used. Fig. ^ illustrates 
the rate of 2 + 1 jet events as function of the invariant mass of the jet system a/S and 
the invariant mass of the whole hadronic final state W. The rate rapidly increases 
towards low V^. Indicated as lines are the phase space boundaries introduced by 
the JADE and the cone algorithm. With a cut on the transverse energy of the jets 
Et > 5 GeV the minimum which can be reached is 10 GeV, since (in the limit of 
massless partons) from = s Zq{l — Zg) follows > Et- When using the JADE 

algorithm with a resolution parameter dij/W"^ ~ s/W'^ ~ ^, low values of \/5 can 
only be reached for low (high) values of W (x). If one takes i/cut = 0.02, \/§ is below 
10 GeV for W < 70 GeV. The JADE algorithm with a high resolution parameter is 
therefore more restrictive than the cone algorithm. In addition to the region of low 
V^, where coUinear and infrared divergences become important, it naturally avoids 
the region of high W (low x) where a large phase space for multi-parton emission is 
available. However, the divergences at small and large Zq can only be avoided by a 
large i/cut- 

Multi-jet production allows a variety of quantitative tests of our understanding 
of QCD dynamics. In a phase space region where the data are well described by 
perturbative QCD, HERA offers the opportunity to extract as over a wide range 
of in a single experiment. This can be done at large x and where the parton 
densities in the proton are well constrained. At low and correspondingly low values 
of ^, gluon initiated processes become increasingly important. Here, the measured 
2 + 1 cross section can be used to determine the gluon density assuming a value 
of as- The sensitivity to both and g{S,-,Q'^) could also be used to simultaneously 
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Fig. 5. a) Fractional transverse energy ^'(r) in a jet defined by a cone algorithm with radius 
R = 1 and 37 < £"( < 45 GeV in the laboratory frame as function of r = Arj"^ + A^^. For 
OPAL jets with Et > 35 GeV were selected, for CDF (DO) the considered Et range of the jets is 
40 (45) < Et < GO (70) GeV. b) ^(r 0.5) as function of the transverse energy of the jet. QGD 
Monte Carlo models are superimposed. They either calculate the dijet cross section in LO and add 
multi-gluon emissions based on the colour dipole model (CDM) or by parton showers (MEPS) or 
treat QCD effects only by parton showers (PYTHIA). 



extract these quantities. Such an analysis is however theoretically and experimentally 
difficult. In the region where the data significantly differ from what is expected by a 
fixed order calculation of perturbative QCD, jet observables are useful tools to gain 
insights in the complex parton dynamics occurring in electron proton collisions. 

2. Jet shapes in DIS 

The internal structure of jets provides useful information on the transition of a 
parton to the complex aggregate of hadrons which can be observed in the detector. 
The capability to describe detailed properties of jets is also an important ingredient 
when comparison of perturbative calculations with experimental data. 

To study the internal structure of jets in the sample of neutral current DIS events 
with Q2 > 100 GeV^ jets are found in the laboratory system using the CDF-cone 
algorithmil with a cone radius R = Ai]'^ + = 1. Only jets are considered which 
have in the laboratory frame a transverse energy Et > 14 GeV and lie in the rapidity 
range —1 < t] < 2. The jet shape is defined as the average fraction of the transverse 
energy of the jets inside an inner cone r concentric to the outer jet cone with radius 
R: 

m = EAr) = I dr'dEr(r')/dr' (4) 



HI preliminary 



Unfolded Data 
MEPJET (A = 600 MeV) 
MEPJET (A = 100 MeV) 
MEPJET (best A) 



10 




0.01 



0.02 



0.03 



0.04 



0.05 



0.06 



0.07 



0.08 



0.09 



0.1 



Fig. 6. Distribution of the resolution parameter j/2 in the JADE algorithm where a 1 + 1 jet event 
switches to a 2 + 1 jet event. Shown arc unfolded data and the next-to-leading order prediction for 
some values of ^qqjj ■ -^dis is the total number of DIS events in the considered kinematic region. 

where Efir) is the transverse energy within an inner cone of radius r and A'jet is the 
total number of selected jets. By definition < \1/ < 1 and '^{R) = 1. The steepness 
of the rise of \l/ describes the collimation of the jet. 

The jet shape \1/ corrected to hadron level for jets with Et > 37 GeV as function 
of r is shown in Fig. 90% of the jet transverse energy is already contained in 
a cone with r = 0.4 around the jet axis. This is very similar to jets produced in 
e~^e~ collisions^ with Et > 35 GeV. The jets selected from pp collisions for 40 (45) < 
Et < 60 (70) GeV measured by CDF (D0)i3ii are considerably broader. This might 
be explained by the fact that in e~^e~ and ep collisions the partons initiating the 
cascade leading to the jet are mainly quarks while in pp collisions they are gluons. 

The jet shape depends on the transverse energy of the jet. Jets get narrower with 
increasing Et as is shown in Fig. 

Monte Carlo models based on a hard subprocess in LO plus additional multi- 
gluon emissions are able to describe the data. Also PYTHIA where QCD effects 
are only treated by parton showers reproduces the main features of the data. The 
best description of the Et dependence of \l/ is obtained when partons are emitted 
according to the colour dipole model. The internal jet structure in DIS at HERA is 
fairly well reproduced by all QCD models. 

3. Determination of the strong coupling constant 

The 2 -|- 1 jet rate defined with the JADE algorithm depends on the jet resolution 
parameter i/cut- For each event y2 is defined as the value of ycut where a 1 -|- 1 jet event 
switches to a 2 -|- 1 jet event. 

A region of phase space is selected where NLO QCD can provide a reasonable 
description of the data, and where effects from multi-gluon emissions are small. In 
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Fig. 7. Left: Values and total error of as{Mz) from various processes. The solid line indicates 
the world average and the band its total error. Right: as{Q) from HERA (open symbols) and 
other processes with increasing Q (closed circles): Fx and oiiad/otot, event shapes and Fhadion/Fiepton 
in e~^e~. The preliminary points using the kr algorithm include the statistical errors only. The 
predictions of the renormalisation group equation for three values of are superimposed as 

lines. 



the HI analysiJi this region is defined by: > 200 GeV^ and > 5000 GeVl 
To avoid events where the invariant mass of the two jet system is small and where the 
cross-section diverges, y2 is restricted to 1/2 > 0.01. The 1/2 distribution, corrected to 
parton level by an unfolding methodic, is shown in Fig. Overlayed is the prediction 
of NLO QCD for three values of as (respectively Aqcd)- In all bins a strong sensitivity 
to as is observed. For one given value of as NLO QCD gives an excellent description 
of the data. The precision which can be expected from this measurement is of the 
order of Aas{Mz)/as{Mz) ^ 5 - 10%. 

Keeping ycut = 0.02 fixed, the number of 2 + 1 jet events (iV2+i) normalized to 
the sum of 1 + 1 and 2 + 1 jet events (A^2+i + ^1+1) can be measured and is sensitive 
to as- Only jets in the acceptance region of 10° < 6jet < 145° and with 2;^ > 0.1 are 
counted. The fractional jet rate R2 = A^2+i/(^i+i + ^^"2+1) can be expressed by: 

R2{Q\ Vent) = A{Q\ y,^t) as{Q^) + B{Q\ y,^t) a^(Q') (5) 

From the coefficients A and B which can be calculated in NLO QCD and the measured 
R21 ocs{Q'^) can be extracted in different bins. Using the renormalisation group 
equation these independent measurements can be used to fit as{Mz). 

For this measurement only the kinematic range of high and high x is considered 
in order to reduce the uncertainties stemming from higher order gluon emissions and 



from the hadronisation of a complex partonic configuration. The theoretically and 
experimentally difficult region towards the proton remnant is excluded. Three values 
of a. extracted by the H10 and ZEUSe3 collaboration are shown in Fig. 0. For 
increasing values of Q^, asiQ"^) decreases as predicted by the renormalisation group 
equation shown for Ajj^ between 100 and 300 MeV (solid lines). An extrapolation 
to as{Mz) yields: 

HI 93 : as{M,) = 0.123 ± 0.012 (stat) ± 0.013 (syst.) 

ZEUS 94 : a,(M^) = 0.117 ± 0.005 (stat) ± 0.005 (exp.) ± 0.007 (th.) 

which is consistent with other values obtained from a large variety of different pro- 
cesses as shown in Fig. ^ (for ref. se The agreement found in e"'"e , ep and pp col- 
lision for reactions at very different scales is an important and successful test of QCD. 
The error on the HERA measurements are even with the present limited statistics 
already competitive. The dominant experimental errors stem from the uncertainty 
on the energy scale for hadrons, the model dependence of assigning partons to jets 
and the phenomenological description of the hadronisation process. Furthermore the 
renormalisation scale and the choice of the input parton density systematically influ- 
ences the result. These results have been obtained with the NLO program PROJET 
or DISJET. 

When particles are combined to jets, the details of how the 4— momenta are added 
are important, in particular, because higher order QCD calculations are only carried 
out for massless partons whereas the jets acquire a mass during the combination 
procedure. In the E-scheme 4-momenta are added to massive jets while in the Eq and 
P-schemes jets remain massless by rescaling either the 3-momenta or the energy (see 
table ID for details). The dependence of the jet cross section on the recombination 
schemes (RS) only appears in the calculation when jets are composed out of more 
than one parton. Corrections of 0{al) due to the different recombination schemes 
are potentially large. When applying different recombination schemes to the data, 
the measured 2+1 jet cross sections defined with the E and Eo algorithms agree 
within few percent, but are lowered by up to 20% for the P-schemeEj. In PROJET 
only the E scheme is implemented. It is therefore not surprising that the extracted 
as{Mz) agrees within 3% when the E and EO scheme is only applied to the data, but 
gives an inconsistent result when using the P schemeS. An additional problem is that 
in the NLO calculation the W algorithm is implemented, while in the experimental 
analysis the JADE algorithm has to be used, i.e. mass terms in dij are neglected. 

The new NLO programs MEPJET and DISENT reveal a strong dependence on 
the recombination schemeil when using the JADE or W algorithm. In addition lead 
these algorithms to large NLO corrections of the dijet cross section. 

The cone or kx jet algorithms seem to be better suited for precision QCD tests, 
since the variation of the cross section when going from leading to next-to-leading 
order is much smaller^l. Moreover, these jet algorithms are less dependent on the 



choice of the renormahsation and factorisation scale and allow a larger phase space 
to be covered. 

The ZEUS collaboration has reanalyzed their data0 using the fcr algorithm. The 
preliminary values of as{Q) obtained in three bins of Q are shown'^in Fig. |^. They 
are consistent with the results obtained with the JADE algorithm. The value extrap- 
olated to the Z° boson mass is: as{Mz) = 0.118 ± 0.008 (stat). 



4. Determination of the gluon density 

The parton distributions are well constrained for larger x where a lot of data 
are available. However, for low x, where gluons play a crucial role, they are not 
well known. In this regime, HERA has the unique possibility to investigate the 
distribution of the gluon density in the proton. The scaling violation of the proton 
structure function Fo offers an indirect way to get a handle on the gluon density. 
This measurementiiii is possible because of the extremely successful description of 
F2 by the NLO DGLAP equations. Observables based on 2 + 1 jet event can - 
by including them in the global analysis - provide a further constraint on the gluon 
density. Dijet rates can be calculated in next-to-leading order and are sensitive to the 
gluon density, since for 10 < g2 < 100 GeV^ the 2 + 1 jet cross section is dominated 
by gluon initiated processes. In this kinematic region less than about 20% of the 2 + 1 
jet events are induced by quarks (see Fig. H). 

An alternative approach is to only base the analysis on the 2 + 1 jet rate. Since 
large invariant jets masses (s > 100 GeV^) are experimentally required to define clean 
jets, the very low x regime can not be reached. These measurements are nevertheless 
important to fill the gap between the results from small-x scaling violations in F2 and 
the data on direct photon and jet production in pp collisions reaching to large-xS. 

The HI collaborationill has already exploited the 2 + 1 jet production to extract 
the gluon density in LO in the range 2 ■ 10~^ ^ ^ ^ 0.18. When processes where only 
one gluon is involved are considered, ^ is directly observable and the determination 
of g{C,,Q'^) is straight forward. This is however only of limited theoretical interest, 
since the LO dijet cross section is subject to large higher order corrections. There 
is moreover a large dependence on the arbitrary choice of the renormahsation and 
factorisation scale. Only, in NLO a reliable prediction of the dijet cross section is 
possible. The 2 + 1 jet cross section, for instance, defined with a cone algorithm for 
jets with Ft > 5 GeV remains constant within 5% when varying the scale by 3 orders 
of magnitude while in LO variations of 40% are calculatedE3. 

A NLO extraction of g{^, Q^) requires, however, an iterative procedure where 
the input gluon density used in the NLO calculation is adjusted until the theo- 
retical dijet cross sections agree with the measurement. This is most efficiently 
done by expanding the gluon density for a discrete set of points {^i,..-,^n) using 

'^Only the statistical errors are included in the figure. 




P, /GeV ' ' ^ , 

Fig. 8. Left: Shape normalised distribution of the transverse energy of jets defined with the 
kt algorithm (in the Breit frame) using /2 as resolution parameter. The data have not been 
corrected for detector effects. Superimposed is a simulation based on LEPTO. Right: The function 
h{x) describing the non-perturbative contribution of the dijet cross section as defined in the HI 
analysis to extract the gluon density in NLO. 



appropriate test functions g{^,Q^) = J2n=T ^i^n,Q'^)Tn{^) such that the integral 
= Yln=T '^(^n, Q^) J / ^ Cqq{x/^, Q"^) Tn{^) in cquatiou D must only be calculated 
once. The coefficient a can then be used in the fitting procedure. Preliminary results 
using an elegant technique based on Mellin momentscd, where the test functions take 
the form T„ = {x/C,)^, have been presented by the HI coUaborationll covering the 
range 0.01 < ^ < 0.1 using the JADE algorithm. 

Another approach also pursued by the HI collaborationlll is to use as test functions 
cubic splineg^d with the property Tj(^j) = 6ij and to introduce the 2 + 1 jet rate as 
additional observable in a global fit analysis minimizing = X^(-^2) + X^iRt+i)- 
This has the advantage that no assumption on the quark densities is needed. The 
correlation of quarks and gluons via the DGLAP equation is naturally taken into 
account. -R2+1 measured in the detector is for each considered kinematic bin 
related to the NLO calculation of the dijet cross section af+i performed by DISENT: 

-^2+1 1 — 7—Tyr\ y^) 



The coefficients MIG2+J correct for the migration of 1 + 1 to 2 + 1 jet events from 
hadron to detector level and can be determined using a fully simulated event sample 
based on LEPTO. The transition from partons to hadrons is modeled by a simple 
ansatz: h^^^'^'^^' [x^Q"^) = (a + /51ogx + 7 log^ x + 51og^ When comparing 
NLO calculations to measurements, usually the data are 'corrected to parton level' 



by assuming that this transition is correctly described in a LO (improved by par- 
ton showers) Monte Carlo simulation program. To introduce a simple hadronisation 
model with additional free parameters in the fit procedure has the advantage that no 
a priori assumption on the underlying parton level is needed to relate the NLO calcu- 
lation to the measurement. Recently, similar models based on power behaved term£i 
have been very successful in describing event shapes in DIsS In this context theoret- 
ical studies suggest that the terms l/Q may also take into account multi-emissions of 
soft giuons associated with the behaviour of the running coupling constant at small 
scalesc3. 

In the kinematic range 

25 < Q2 < 100 GeV^ and 2 ■ 10"^ < x < 2 • lO^^ jg^g 
are selected using the kx algorithm with dij = 0.5 and as scale. The jet rate for 
jets with 7° < 9jet < 110° and Zg > 0.15 is measured. LEPTO is able to describe the 
shape of most distributions associated to the hard subprocess. As an example Fig. |] 
shows the Et distribution of the jets as measured in the detector. LEPTO, however, 
fails to describe the absolute jet rate. The measurement is always about 20% above 
the predicted jet rate. This is true in the whole considered {x,Q'^) plane, but most 
pronounced for low (10 < < 14 GeV^) and low x (2 ■ 10"^ < x < 1 ■ lO'^). 
The observed difference can not be accounted for by systematic effects. The most 
important systematic error of about 3 — 13% on the dijet rate is mainly due to the 
uncertainty of the hadronic energy scale. 

The fit to the F2 data and the jet rates leads to a good x^(-R2+i) = 10 for 11 
degrees of freedom, but no influence of the jet data on the extracted parton density 
functions is found in this analysis. Also the error band on the extracted gluon density 
is only marginally improved. This means that the extracted gluon density is essen- 
tially constrained by the inclusive F2 data. The fitted non-perturbative function h{x) 
is shown in Fig. |^. Without a strong x dependence of h{x) no satisfactory fit can be 
performed. 

A possible explanation of this results is that the way in which the kx jet algorithm 
is used in this analysis leads to jets with relatively low transverse energy of 3 — 
5 GeV (see Fig. H) such that the found jet system does not always stem from the 
hard subprocess. Another explanation is that a fixed order calculation is not sufficient 
to account for the data in this phase space region. The difference between the NLO 
calculation and the data has therefore to be absorbed in the non-perturbative function 
h{x) and consequently no sensitivity of the jet rates to the gluon density is found. 



5. Comparision of jet rates to NLO-QCD 

A similar observation about the incapability of NLO QCD to describe the dijet 
rates has been made by the ZEUS collaborationEl. In this study jets are selected using 
a cone algorithm in the laboratory frame requiring > 4 GeV in the hadronic center 
of mass and the laboratory frame and excluding the forward region by r^iab < 2. The 
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Fig. 9. 2 + 1 jet rate cross section in function of the momentum fraction of the parton entering the 
hard subprocess ^. Data have been corrected to parton level using the QCD Monte Carlo program 
LEPTO 6.3. Superimposed is a prediction from the QCD matrix element in LO and in NLO in as- 



data are then corrected to parton level using LEPTO 6.3. The correction factors are 
between 1.2 and 1.8 for 4 ■ lO^^ < ^ < 8 ■ 10"^ but get very large (up to 3 — 4) 
for ^ around 10^^. The shape of variables connected to the hard subprocess like 
Et, Tj and ^ are well described by NLO programs. The absolute value of the 2 + 1 
cross section is, however, found to be about 30% higher in the data. This large 
difference cannot be explained by experimental effects like variation of the energy 
scale or the model dependence in the correction procedure nor by varying theoretical 
choices like different parameterizations of the parton densities or different factorisation 
and renormalisation scales. 

That NLO QCD is not able to describe the dijet rates has recently been reported 
by the HI collaboration@0. The dijet rate is measured with a cone algorithm using a 
cone radius R = 1 and requiring Et > 5 GeV in the hadronic center of mass system. 
The rapidity difference between the two jets has to be Ar/ < 2. The 2 + 1 jet rate 
corrected for detector effects is shown in Fig. 10 as function of x and Q^. NLO QCD 
(DISENT) as well as a conventional QCD model based on the exact LO coefficient 
functions plus additional DGLAP parton shower (RAPGAP 'DIR') clearly fail to 
describe the data. This is true for the whole kinematic range 5<Q^ < 100 GeV^ and 
lO"'^ < X < 10~^, but is most pronounced in the low x and region. In this region 
the prediction of the NLO QCD calculation is a factor of 2 higher than the QCD 
model, but is by the same factor below the data. Varying the input parton density 
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Fig. 10. 2 + 1 jet rates as function of and x corrected to hadron level. Jets are defined by 
the cone algorithm with Et > 5 GeV in the hadronic center of mass system. The difference in 
pseudorapidity between the two jets has to be smaller than 2. The band indicates the systematic 
error. Superimposed are the predictions from ARIADNE, a next-to-leading calculation (DISENT) 
and from RAPGAP (the direct and resolved contribution are shown separately). 



functions leads to a change in the jet rate by up to 30%, but cannot account for the 
large discrepancy with the data. 

The CDM model as implemented in ARIADNE is able to describe the data. In 
this model gluons are emitted from a chain of colour dipoles spanned between the 
current quark and the proton remnant. The first emission in the cascade is corrected 
by the matrix element to first ordeiS. It has been observediil that quark induced 
processes have a up to a factor of 8 higher cross section than expected from conven- 
tional QCD, because the parton distribution functions in equation |I] are replaced by a 
phenomenological expression better adapted to the dipole picture. Partons with short 
wavelengths having high transverse momenta (Pt) only resolving part of the extended 
colour charge distribution of the proton remnant are suppressed by (/i/Pt)" where 
/i and a are phenomenological parameters usually set to 1. A similar suppression 
depending on the virtuality is introduced in the direction of the scattered quark. 
It has been argued that the partonic state as encountered in ARIADNE is closely re- 
lated to the one expected from a BFKL scenarioil. In ARIADNE the hardest parton 
emission is not bound to the photon vertex, but can happen anywhere between the 
photon and the proton remnant. 

Such a feature can also be artificially introduced in the conventional QCD model 
by assuming that the photon emitted from the incoming electron can fluctuate into 
quarks and gluons before entering the hard subprocess. This can be modeled by 
assuming a structure in the virtual photon like e.g. parameterized in the SaS modeS. 
Fig. p!0| shows that adding this contribution (RAPGAP 'RES') to the processes where 
the photon is assumed to be point-like (RAPGAP 'DIR') does also reproduce the 
measured jet rate. 

Only models where the hardest parton emission is not bound to occur at the 
photon vertex are able to describe the data. This indicates that in the phase space 
region tested by 2-1-1 jet production higher order processes which cannot be described 
by NLO QCD are important. Whether it is sufficient to simply add processes in the 



photon direction without being in conflict with other HERA data, is not yet clear. 
The phenomenological treatment of higher orders in the colour dipole model has been 
proven to be extremely successful in describing all the data on the hadronic final state 
which have been measured at HERA^'i. It remains however a puzzle why the shape 
of all distributions associated with the hard dijet system, like Et, Xp, t] etc., is well 
described by conventional QCD. It is hard to imagine that the importance of higher 
orders only affects the normalisation and not the shape of the distributions. 

6. Conclusions 

The understanding of jet production in DIS remains one of the main physics goals 
at HERA. While with the first HERA data, topics like the measurement of the strong 
coupling constant and the gluon density have been attacked with great enthusiasm, 
it turned out that these goals can only be achieved in a (very) restricted phase space 
region. 

Existing data on jet production have always been tailored to a specific goal and 
present themselves today as a loose collection of scattered information. While con- 
ventional QCD is able to describe the data for > 100 GeV^ and 6j(,t > 10° using 
the JADE algorithm, a discrepancy of about 30% is found for the low region 
when using the cone or the /ct algorithm and excluding the forward region. The fully 
acceptance corrected jet rates in this region disagree by about a factor of two. 

It seems that before precisely pinning down the strong coupling constant and the 
gluon density, one has to systematically investigate the applicability of NLO QCD in 
terms of kinematic phase space and the way how the jet observables are defined (jet 
algorithm, hard scale etc.). 

In future, it will be necessary to measure jet cross sections corrected for detector 
effects extending to very low and very large x and in order to find out where 
perturbative QCD is able to describe the data and where the mechanism governing 
jet production are more complicated. To achieve this goal the understanding of QCD 
models will play a major role. 

7. Acknowledgements 

I would like to thank my colleagues J. Dainton, G. lacobucci, M. Kuhlen, E. Mirkes 
and M. Wobisch for their critical reading of the manuscript and helpful comments. 

8. References 

1. A. Mendez, Nucl. Phys. B145 (1978) 199; R. Peccei and R. Riickl, Nucl. Phys. B162 (1980) 
125; C. Rumpf, G. Kramer and J. Willrodt, Z. Phys. C7 (1981) 337. 

2. Y.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641; V.N. Gribov and L.N. Lipatov, Sov. J. 
Nucl. Phys. 15 (1972) 438 and 675; G. Altarelli and G. Parisi, Nucl. Phys. 126 (1977) 297. 



3. A. Edin, G. Ingelman and J. Rathsman, Z. Phys. C75 (1997) 57. 

4. K. Flamm, Phd thesis, University of Hamburg, DESY FHl-96-04 Hamburg, 1996. 

5. D. Graudenz, Phys. Lett. B256 (1992) 518; Phys. Rev. D49 (1994) 3291; Comp. Phys. 
Comm. 92 (1995) 65. 

6. T. Brodkorb and J. G. Korner, Z. Phys. C54 (1992) 519; T. Brodkorb and E. Mirkes, Z. 
Phys. C66 (1995) 141. 

7. JADE CoUab., W. Bartel et aL, Z. Phys. C33 (1986) 23. 

8. E. Mhkes and D. Zeppenfeld, Phys. Lett. B380 (1996) 205. 

9. K. Fabricius, G. Kramer, G. Schierholz and I. Schmitt, Z. Phys. Cll (1981) 315; G. Kramer 
and B. Lampe, Fortschr. Phys. 37 (1989) 161; W.T. Giele, E.W.N. Glover, Phys. Rev. D46 
(1992) 1980: W.T. Giele, E.W.N. Glover and D. A. Kosower, Nucl. Phys. B403 (1993) 663. 

10. S. Cataui and M. Seymour, Nucl. Phys. B485 (1997) 291. 

11. R.K. Ellis, D.A. Ross and A.E. Terrano, Nuel. Phys. B178 (1981) 421. 

12. S. Catani and M. Seymour, Phys. Lett. B378 (1996) 287. 

13. S. Catani, Y.L. Dokshitzer and B.R. Webber, Phys. Lett. B285 (1992) 291. 

14. E. Mirkes and D. Zeppenfeld, in: 5th Int. Workshop on DIS and QCD, eds. J. Repond and 
D. Krakauer, Chicago (USA) 1997. 

15. L. Lonnblad, Comp. Phys. Comm. 71 (1992) 15. 

16. G. Ingelman, in: Physics at HERA, eds. W. Buchmiiller and G. Ingelman, Hamburg (Ger- 
many) 1991, VoL 3, p. 1366; G. Ingelman, A. Edin and J. Rathsman, Comp. Phys. Comm. 
101 (1997) 108. 

17. G. Marchesini et al., Comp. Phys. Comm. 67 (1992) 465. 

18. T. Sjorstrand, Comp. Phys. Comm. 82 (1994) 74. 

19. H. Jung, Comp. Phys. Comm. 86 (1995) 147. 

20. M. Bengtson and T. Sjostrand, Z. Phys. C37 (1988) 465; M. Bengtson, G. Ingelman, T. 
Sjostrand, in: HERA Workshop 1987, ed. R.D. Peccei, Hamburg (Germany) 1988, Vol. 1, 
p. 149. 

21. G. Gustafson and U. Petterson, Nucl. Phys. B306 (1988); G. Gustafson, Phys. Lett. B175 
(1986) 453; B. Andersson, G. Gustafson, L. Lonnblad and U. Petterson, Z. Phys. C43 (1989) 
625. 

22. B.R. Webber, Nucl. Phys. B238 (1984) 492. 

23. B. Andersson et al., Phys. Rep. 97 (1983) 31. 

24. T. Sjostrand, Comp. Phys. Comm. 39 (1986) 347; T. Sjostrand and M. Bengtsson, Comp. 
Phys. Comm. 43 (1987) 367, and for JETSET 7.3, T. Sjostrand, CERN-TH-6488-92 (1992). 

25. CDF Collab., F. Abe et al., Phys. Rev. D45 (1992) 1148. 

26. OPAL Collab., R. Askers et al., Z. Phys. C63 (1994) 197. 

27. CDF Collab., F. Abe et al., Phys. Rev. Lett. 70 (1993) 713. 

28. DO Collab., S. Abachi et al., Phys. Lett. B357 (1995) 500. 

29. M. Weber, in: 5th Int. Workshop on DIS and QCD, eds. J. Repond and D. Krakauer, 
Chicago (USA) 1997. 

30. V. Blobel, DESY 84-118, Hamburg 1984, and CERN Comp. School 1984. 

31. HI Collab., T. Ahmed et al., Phys. Lett. B346 (1995) 415. 

32. ZEUS Collab., M. Derrick et al., Phys. Lett. B363 (1995) 201. 

33. Particle Data Group, R.M. Barnett et al., Phys. Rev. D54 (1996). 

34. K. Rosenbauer, in: Workshop on DIS and Related Phenomena, ed. G. D'Agostini and A. 

Nigro, Rome (Italy) 1996, p. 444. 

35. E. Mirkes and D. Zeppenfeld, in: Workshop on DIS and Related Phenomena, ed. G. 
D'Agostini and A. Nigro, Rome (Italy) 1996, p. 428. 

36. E. Mirkes and D. Zeppenfeld, Acta Phys. Polon. B27 (1996) 1393. 

37. T. Trefzger, in: Workshop on DIS and Related Phenomena, ed. G. D'Agostini and A. Nigro, 
Rome (Italy) 1996, p. 434. 



38. HI Collab., S. Aid et al., Nucl. Phys. B470 (1996) 3. 

39. ZEUS Collab., M. Derrick ct al., Phys. Lett. B345 (1995) 576. 

40. R.C.E. Dcvcnish and A.D. Martin, in: Workshop on DIS and Related Phenomena, ed. G. 
D'Agostini and A. Nigro, Rome (Italy) 1996, p. 193. 

41. HI Collab., S. Aid et al., Nucl. Phys. B449 (1995) 3. 

42. D. Graudenz et al. Z. Phys. C70 (1996) 77. 

43. F. Zomer, in: 5th Int. Workshop on DIS and QCD, eds. J. Repond and D. Krakauer, 

Chicago (USA) 1997. 

44. G. Lobo, in: 'Future Physics at HERA', eds. G. Ingelmann, A. De Roeck and R. Klanner, 
Hamburg (Germany) 1996, p. 537. 

45. Y.L. Dokshizer and B.R. Webber, Phys. Lett. B352 (1995); R. Akhoury and V.L Zakharov, 
Phys.Lett. B357 (1995) 646. 

46. HI Collab., C. Adloff et al., DESY 97-098, Hamburg, 1997. 

47. D. Mikunas, in: 5th Int. Workshop on DIS and QCD, eds. J. Repond and D. Krakauer, 
Chicago (USA) 1997. 

48. J. Spickermann, in: XXXIInd Rencontres de Moriond, ed. J. Tran Thanh Van, Les Arcs 
(France) 1997. 

49. M. Wobisch, in: 5th Int. Workshop on DIS and QCD, eds. J. Repond and D. Krakauer, 
Chicago (USA) 1997. 

50. L. Lonnblad, Z. Phys. C65 (1995) 285. 

51. J. Rathsman, Phys. Lett. B393 (1997) 181. 

52. L. Lonnblad, Z. Phys. C65 (1995) 285; A. H. Mueller, Nucl. Phys. B415 (1994) 373; J. 
Rathsman, Phys. Lett. B393 (1997) 181. 

53. G. Schuler and T. Sjostrand, Z. Phys. C68 (1995) 607 and Phys. Lett. B376 (1996) 193. 

54. N. Brook et al., in: Future Physics at HERA, eds. G. Ingelman, A. DeRoeck and R. Klanner, 
Hamburg (Germany) 1996, Vol. 2, p. 613. 

55. T. Carli, in: Workshop on Deep-Inelastic Scattering and Related Phenomena, ed. G. 
D'Agostini and A. Nigro, Rome (Italy) 1996 p.415. 



